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Abstract. By means of quantum stochastic calculus we construct a model for 
an atom with two degenerate levels and stimulated by a laser and we compute 
its fluorescence spectrum; let us stress that, once the model for the unitary 
atom-field dynamics has been given, then the spectrum is computed without 
further approximations. If only the absorption/emission term is included in the 
interaction, we reobtain the Mollow spectrum in the case of a monochromatic laser 
and the Kimble-Mandel spectrum in the case of a "phase diffusion model" for a 
non monochromatic laser. However, our model can describe also another type of 
light scattering, a "direct scattering" due to the response of the atom as a whole, 
which we expect to be small, but which interferes with the scattering due to the 
absorption/emission channel. When both the scattering channels are introduced 
we obtain a modification of the Mollow-Kimble-Mandcl spectrum, which shares 
the main features with the usual case, but which presents some asymmetries even 
in the case of no detuning. 



PACS numbers: 42.50.-p, 42.50.Ct, 42.50.Hz, 02.50.Ey 

1. Introduction 

Quantum stochastic calculus (QSC) jjj— Q is a noncommutative analogue of the 
classical Ito's calculus and it has had many applications in modelling quantum 
optical systems Q-Q. In particular, inside QSC a Schrodinger equation for unitary 
operators can be formulated, the Hudson-Parthasarathy equation (|5]) jl| , which is well 
suited to treat the interaction between matter and electromagnetic field within the 
usual approximation schemes of quantum optics p"^ ] . 

The simplest, non trivial matter-field system is a two-level atom stimulated by 
a laser. The fluorescence spectrum of such a system, in the case of a perfectly 
monochromatic laser, was obtained by Mollow pl| ; then, the case of a laser with 
a Lorentzian spectrum was developed by Kimble and Mandel jl6|, [l7|. QSC and 
Hudson-Parthasarathy equation give not only the way of giving an unified treatment 
of such a system and of connecting it to the photon-detection theory, but they allow 
also to explore possible corrections to the dynamics and to obtain modifications of the 
Mollow spectrum. In the usual treatments, the scattering of the light by the atom 
is described by the absorption/emission process due to the electric dipole interaction 
involving only two states of the atom. But even if the absorption/emission process 
would be forbidden and the atom would be frozen in the up or down level, some 
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scattering of light would remain, due to the response of the atom as a whole; we can 
call it "direct scattering" . For instance, in a perturbative development in Feynman 
graphs, scattering processes would be generated also by virtual transitions starting and 
ending in one of the two states left in the description, but involving as intermediate 
states the other ones, which have been eliminated in the final description. The use 
of Hudson-Parthasarathy equation allows, even when the atom is approximated by 
a two- level system, for the introduction of both the "absorption/emission channel" 
and the "direct scattering channel" . This modified model, in the case of a perfectly 
monochromatic laser, has been presented in and developed in |l{|; the first results 
for a non monochromatic laser have been obtained in |2C[| . 

In Section || we give the mathematical model for the dynamics of the two-level 
atom interacting with the electromagnetic field, including the polarization of the light 
and the degeneracy of the up and down levels, which has never been considered before 
and which has an effect on the angular dependence of the spectrum. We consider the 
usual dipole interaction and the possible corrections giving the direct scattering; the 
so called phase diffusion model is used for the laser as in . In Section || the general 
formulae for the spectrum are obtained and its angular dependence is discussed. Then, 
the total spectrum, integrated over the whole solid angle, is studied. In the case of 
only the absorption/emission process in the interaction term, we reobtain the Mollow 



spectrum [15 for a perfectly monochromatic laser and the Kimblc-Mandel spectrum 
Jl6| for a non monochromatic laser with Lorentzian spectrum. In the case of the 
modified model, when also the direct scattering is included, but the laser is perfectly 
monochromatic, we reobtain the results of Q, although there the level degeneracy 
was not taken into account. According to the values of the various parameters, the 
fluorescence spectrum can present a well resolved triplet structure as in the Mollow 
spectrum, but also it can be distorted and made asymmetric by the presence of the new 
terms. Experiments pl[~|25[ confirm essentially the triplet structure; some asymmetry 
has been found and its origin has been attributed to various causes. The presence of 
the new terms we have introduced could give a model for treating such an asymmetry. 
When the laser bandwidth is taken into account, the modified model gives corrections 
to the spectrum found in |l6| by Kimble and Mandel. 

2. The model 

2.1. The unitary atom-field dynamics 

2.1.1. Fock space and the Hudson-Parthasarathy equation. Let us recall the rules of 
QSC and the Hudson-Parthasarathy equation; this is just to fix our notations, while 
for the proper mathematical presentation we refer to the book by Parthasarathy Q . 
We denote by T the Boson Fock space over the "one-particle space" Z®L 2 (R + ) , 



where Z is a separable complex Hilbert space; we shall see in section |2.1.3| how to 
choose Z. A vector / in Z <E> L 2 (M + ) ~ L 2 (M. + ;Z) is a square integrable function 
from R + into Z. Let {hk, k > 1} be a c.o.n.s. in Z and let us denote by Ak(t), 
Af.(t), Aki(t) the annihilation, creation and gauge (or number) processes associated 
with such a c.o.n.s.; we can write formally Ak(t) = J*ak(s)ds, A\,(t) = f*a\,(s)ds, 
Afcj(t) = fl a k (s)ai(s) ds, where ak(t) 7 a\{t) are usual Bose fields, satisfying the 
canonical commutation rules [a/c(i), a\{s)\ = SkiS(t — s), . . .. The Fock space T is 
spanned by the coherent vectors e(/), / 6 L 2 (R+; Z), whose components in the 
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0, 1, . . . , n, . . . particle spaces are 

e(/) = exp HH/H 2 ) (l, /, (2!)- 1 / 2 / »/,..., (d)" 1 / 2 /®", . . .) ; (1) 

in particular the vector e(0) is the Fock vacuum. The coherent vectors ([!]) are 
normalized and they are the eigenstates of the annihilation operators; indeed, 
we have A k (t)e(f) = J*(h k \f(s))dse(f). In QSC integrals of "Ito type" with 
respect to dAk(t), dA\.(t), dA k i(t) are defined. The main practical rules to 
manipulate "Ito differentials" are the facts that dA k {t), dA' k (t), dA k i(t) commute with 
anything containing the fields only up to time t (objects called adapted processes) 
and that the products of the fundamental differentials satisfy dA k (t) dAj (t) = 
5 k idt, dA k (t) dA rl (t) = SkrdMt), dA kr (t)dA](t) = S H dAl(t), dA kr (t) dA s i(t) = 
5 rs dAki(t); all the other products and the products involving dt vanish. 

Let Hbea separable complex Hilbert space (the system space) and let R k , k > I, 
Skh k,l > 1, H be bounded operators in H such that W = H, J2 k R k Rk 1S strongly 
convergent to a bounded operator, and 

^S k i®\h k ){h\=:SeU{H®Z), (2) 

kl 

U(*) "means unitary operators on *"; we set also 

R:H^H<E)Z, R£ := ^(Rkt) ® h k , (3) 

fc 

K = H- 1 -Y,R{Rk = H- 1 -R^R. (4) 

fc 

Then (see reference ^ Theorem 27.8 p. 228) there exists a unique unitary operator- 
valued adapted process U(t) satisfying {7(0) = II and 

dC/ W = Rk dA lW + J2 (Ski - 5ki) dA kl {t) - R{Sm dMt) - iKdt^U(t).(5) 

The Bose fields introduced here represent a good approximation of the electro- 
magnetic field in the so called quasi-monochromatic paraxial approximation [ p6| p"2| . 
Now, T is interpreted as the Hilbert space of the electromagnetic field; A\ (t) creates 
a photon with state h k in the time interval [0,t], A k (t) annihilates it, A k k(t) is the 
selfadjoint operator representing the number of photons with state h k in the time 
interval [0, £]. The operator U(t) represents the evolution operator for a compound 
system (say an atom plus the electromagnetic field), in the interaction picture with 
respect to the free dynamics of the field. 



2.1.2. The two-level atom. The experiments 21 - |25[ involve or the hyperfinc 
component of the D2 line of sodium with up level 2P 3 / 2 , F = 3 and down level 
25i /2) F = 2, or the levels 6s6p ^1, F = 1, and 6s 2 S , F = of 138 Ba; F is the 
total angular momentum. 

In order to describe an atom with two degenerate levels as in the experimental 
situation, we take 

H = C 2F - +1 ©C 2F++1 , F+=F_ + 1, (6) 

where F_ is integer or semi-integer. We denote by \F± , Af), M = —F±,...,F±, 
the angular momentum basis in 7i\ the parities of the states of the two levels 
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must be opposed, let us say they are e±, with e + e_ = — 1. Let us denote by 
P± = Yl^a=-F± ' M)(F± " M\ the two projection on the up or down states. Then 
the free energy H, contained in the quantity (|j), must be given by 

H=±lo (P + -P_), ^ >0, (7) 
where the atomic frequency u>o must already include the Lamb shift. 

2.1.3. The photon space. In the approximations we are considering |2^, the 
fields behave as monodimensional waves, so that a change of position is equivalent to 
a change of time and viceversa. Then, the space Z has to contain only the degrees of 
freedom linked to the direction of propagation and to the polarization. To describe a 
spin-1 0-mass particle we use the conventions of Messiah |^7j pp. 550, 1032-1037. 

A c.o.n.s. in Z is given by \j,m;vj) = Qf m , j = 1, 2, . . ., m = -j, -j + 1, . . . , j, 
w = ±1; jm is the total angular momentum, w — +1 denotes electrical multipole, 
zu = — 1 denotes magnetic multipole and (— ly-uj is the parity. By using the spherical 
harmonics Yi m (9, <f>) and the orbital angular momentum operator £, one has 

ej„\ = J— £Y jm , = ip x 9Ti , (8) 

where p is the direction versor given by p\ = sin# cos</>, p2 = sin 9 sm.<f>, p3 = cos 9. 

2.1.4- The interaction. Let us consider now the terms with the creation and 
annihilation operators in the dynamical equation (||); they must describe the 
absorption/emission process. By asking spherical symmetry, parity conservation and 
only electrical dipole contribution in R, we must have 

R = a \e+;F-;l;F + ,M)(F + ,M\, aeC, a^O, (9) 

M=-F + 
F- 1 

|e+;F_;l;F + ,M) = ^ £ |F_, mi) ® |1, m 2 ; +1)(F_, m i; 1, m 2 \F + , M), (10) 

mx=—F- m2 — — l 

where by (F-, mi; 1, rri2\F + , M) we denote the Clebsch-Gordan coefficients ( |p7|| pp. 
560-563). The interaction terms containing R are responsible for the spontaneous 
decay of the atom and, as we shall see, \a\ 2 turns out to be the natural line width. 
Note that 

R k =P-R k P+, R^R=\a\ 2 P+. (11) 

The interaction term containing the A-process must give the residual scattering 
when the atom is frozen in the up or down level, so we take the unitary operator S of 
the form 

S = P+ <g> S+ + P_ O S~ , S ± e U{Z) . (12) 
Then, by spherical symmetry and parity conservation we must have 

oo j 

s ± = ^2 X! X! b'' 771 ;^) 6 ^! 2 ^^;^)}^^;^ 7 !: ( 13 ) 

j — l m——j tn— ±1 

< 6±(j;ru) < 2tt; the <5's are the phase shifts for the atom frozen in the up or 
down level. Quantities like loq, a, S±(j;w) are phenomenological parameters, or, 
better, they have to be computed from some more fundamental theory, such as some 
approximation to quantum electrodynamics. 
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2. 2. Some dynamical features of the model 

2.2.1. The master equation. Equation (|^) already contains the usual Markov 
approximations; indeed, when the initial state of the field is a coherent vector, then the 
atomic reduced dynamics turns out to be given by a master equation without need of 
introducing further approximations jl], |], [l^] . If po is the initial atomic density matrix 
on Ti and the coherent vector e(/) is the state of the field, then the atomic reduced 
statistical operator 

p(t;f) = Tr, {U(t)( Pa ® |e(/))<e(/)|)E/(i)t} (14) 

turns out to satisfy the master equation 

^ p(t; /) = r(/(t)) [p(t;f)], (15) 
where the time dependent Liouvillian turns out to be 
£(/(i))[p] = -i [H(f(t)),p] 

+ lY,{[ Rk (f( t ))p>Mf( t )y] + [MW))>pMf (*))*])> ( 16 «) 

k 

Rk(f(t)) =R k + (h k \S+f(t))P+ + (h k \S-f(t))P- , (166) 
H(f(t)) =H + l -Y J {{S-f{t)\h k )R k - (h k \S-f(t))Rl). (16c) 

k 

Note that, when the field is in the Fock vacuum, only spontaneous emission is present; 
indeed, for f(t) = 0, the Liouvillian ( |l6o| ) reduces to 

i 1 1 

£(0)[p}=- 1 - lo [P + -P^,p]--\a\ 2 (P + p + pP + ) + \a\ 2 £ A m pA] n , (17a) 

m — — 1 

A m = 2J toi)(F_, mi; 1,771^+, mi +m)(F + , ?7ii +m| , (176) 

mi=-F_ 

which describes decay according to the usual selection rules for electric dipole and 
fixes the meaning of \a\ 2 . 



2.2.2. The balance equation for the number of photons. Let us introduce the 
observables "total number of photons in the time interval [0,7;]" before and after the 
interaction with the atom 

7V in (t) = ^A fcfc (i), N° ut (t) = tf(t)N in (t)U(t). (18) 

k 

By using the computations on the increments of the "output fields" given in fl^| , we 
obtain in the present model the balance equation 

N™(t) + \{P+- P-) = N° ut (t) + 1 U(tf(P + - P_)U(t) . (19) 

Such an equation says that the number of photons entering the system up to time 
t plus the photons stored in the atom at time is equal to the number of photons 
leaving the system up to time t plus the photons stored in the atom at time t. An 
equation for mean values similar to (|l^) has been used for restricting the possible 
forms of the interaction term in the simpler model studied in [[19[ . 
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2.3. The laser and the equilibrium state 

2.3.1. The phase diffusion model for the laser. As in |]l6| , [TtJ we take as the state of 
the laser a mixture of coherent vectors E[|e(/))(e(/)|] with 

f{t) =expj-i(wt + VBW{t)y\e{T -t)X, XeZ, lu > , B>0, (20) 

W(t) is a standard Wiener process, 9 is the usual step function and T is a large 
time; T — * +oo in the final results is always understood and E means the classical 
expectation over the Wiener process. Note that this model of the laser implies a 
Lorentz spectrum of bandwidth B; indeed, for T — > oo, one obtains 



2nJoo ~ E[(/(t)|/( f + r))]dr = H|Ar (|/ _^ + ^ /4 . (21) 

The whole model is meaningful only for u> not too "far" from w - 

In order to describe a well collimated laser beam propagating along the direction 
z (6 = 0) and with right circular polarization, we take 

A = afl c iS X+ , Q > , S e [0, 2tt) (22) 

U0, <t>) = , 1[0|A ' l(g) (~) (i + ij) , (23) 

A0^37r(l-cosA0) \ V2J\ J J' V ' 

where l[o,Ae](^) = 1 for < < A9, l[o,Ae](^) = elsewhere; in all the physical 
quantities the limit A9 J, will be taken. Note that the power of the laser 
M|A|| 2 = §Ma| 2 f27(A60 2 diverges for A9 J. 0, because we need a not vanishing 
atom- field interaction in the limit. In the following we shall need the relation 



0', m; zu\X + ) = -mJ ^ ^ 1 8 m . 1 (24) 

The choice of a circular polarization is done in some experiments, because in this 
way, in the long run, only two states are involved in the dynamics ( |fj3f p. 206) and the 
usual theory has been developed for two-states systems [O, |l6) ; we shall see that also 
in our modified model the right circular polarization pushes the atom in the states 
\F_,F_),\F + ,F + ). 



2.3.2. The equilibrium state. By a stochastic unitary transformation, we obtain from 
(Q) and (|2^) an atomic density matrix 

p(t) = exp (ut + VBW(tj) (P+ -P-)\ p(t; f) 

x exp{-i(W\/W(i))OP+- J P-)}, (25) 
which satisfies the "Ito's" master equation 

dp(t) = C B [p(t)]dt + l - VB[P+ - P- 1P {t)]dW(t) (26) 
and whose mean value satisfies the ordinary master equation 

p(t)=E[p(t)], ^ t p(t)=£ B [p(t)]. (27) 
The time independent Liouvillian entering equations ( [26] ) and ( p7| ) is given by 
L B [p] = C [p] + | ( [(P+ -P-)p,P+-P-] + [P + -P-,p(P+-P-)]), (28a) 
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2 

k 



H = \{u Q - u)(P+ -P_) + lJ2 (((S- + l)X\h k )R k - (h k \(S- + 1)A)4 



k 



+ i «A|S+A) - (5-+A|A>) P + + i ((A|5-A) - (5"A|A)) P_ , (28c) 
where -D(/j) is a kind of effective dipole operator 

D{h) = J2(h\h k )Rk + (h\ {S+ - il) X)P+ + (h\ (S- - 11) A)P_ 

l 

= a ^ (h|l,m;+l)A m + afte i ^(/i|(S e -fl)A + ).P £ . (29) 

m— — 1 e— ± 

Now, let us set 

P x = 1- |1,1;+1)(1,1;+1| , e = Im(5 + A + |Pj_S-A + ) , (30a) 

,9± = Pl (S ± - H)A+ , Ag = , 0+ - g _ , (306) 

8+ = 8+(l;+l), s = 8+ -6-, (30c) 

We have to ask |e| < +oo, \\g±\\ < +oo, for A6 J. 0; roughly speaking, 5± must 
introduce small corrections even when the norm of A+ diverges. Then, we have 

P+C B [p\P+ = -\a\ 2 P +P P+ - L^P-pP+ - P+pP-L- , (31a) 
l 

P-C B [p]P- = \a\ 2 AnP + pP+A} n + L-p + pP_ + P-pP+L-^ , (316) 

m=— 1 

P-£ B [p]P+ = (P+C B [p f ] P-Y = L+P +P P+ - P_pP_L- 

+ [i(uj Q -LU + Im{S + X\S-\)) - i (H 2 + B + ||(5+ - S-)A|| 2 )] P-pP+, (31c) 

L± = J2^ S±X \ h ^ Rk = -\ lapfie-^+^Ai , (32a) 

k 

(5* - 11) A+ = -ie i,5 ±sin5 ± |l,l;+l)+5 ± , (326) 

Im(5+A|5-A) = |a| 2 ^ 2 (e - i sin 2s) , (32c) 

|| (S+ - S-) Af = \a\ 2 n 2 (sin 2 s + || A ff || 2 ) , (32d) 

Then, one can check that the equilibrium state poo for the master equation ( p7|) 
is supported by the two extreme states 

\1) = \F+,F+), |2) = |F_,F_>. (33) 
By using these two states the density matrix is given by 

/ Q 2 d 1 nexp[i(8 + 2S^)]d 2 

Poo = | , (34) 

\ exp [-i (<J + 2<5_ )] d 3 1 - Q 2 di 

*\ 1 / 

d = ( d 2 , Gd = to , w = - 1 | , (35) 

* / V 1 
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-1/2 -1/2 \ 



G = 



Sl 2 e ls cos s 



y fl 2 e 1S cos t 



b J 



b = \ [1 + y + O 2 (|| A<?|| 2 + sin 2 s)] - i (z + \n 2 sin 2s) , 



z = z — Q 2 e , z = (cj — cJq) /(o^ 



V = B/\a\ 



It turns out 



d\ = Re d , d2 = d, d^ = d, 

1 + y + ti 2 (||A5|| 2 + sin 2 s) +i(2z-ft 2 sins coss) 
= 4P + T 2 



r 2 = + y + Q 2 || A ff || 2 ) 2 + fl 2 (2 + 2y + 2ft 2 || A ff || 2 + ft 2 si 



sm 2 s 



(36) 

(37) 
(38) 

(39) 
(40) 

(41) 



3. The fluorescence spectrum 

3.1. Heterodyne detection 

3.1.1. The output current. A way to obtain the spectrum of the stimulated atom is by 
means of balanced heterodyne detection] in such a scheme the light emitted by the atom 
is made to interfere with the light of a strong monochromatic laser (the local oscillator) 
and then detected by a couple of photocounters |2§] . A mathematical description of 
such a detection scheme has been given inside the theory of continuous measurements 
p2l pjl ; the output current of the detector is represented by the operator 



I{y,h;t) 



F(t-s)e ivs dA(h;s) + he. 



dA(h;t) 



E 

k 



where F(t) is the detector response function, say 



F(t)=k 1 i 



exp 



7 



7>0, 



(h\h k )dA k (t) , (42) 



(43) 



fci 7^ has the dimensions of a current, v is the frequency of the local oscillator and 
h £ Z, \\h\\ = 1. Any information on the localization and on the polarization of the 
detector is contained in h. We assume that the detector spans a small solid angle, so 
that h is given by 

1at(0W) ^ ^ 



(44) 



^AT| 

where AT is a small solid angle around the direction (0,0), AT | {(9,<fr)}, |AT| ~ 
sm8dd d<f> and e is a complex polarization vector, \e(9' , <f/)\ — 1. Moreover, we assume 
that the transmitted wave does not reach the detector, i.e. 6 > and so 



(h\X) =0. 



(45) 



From the canonical commutation relations for the fields one has that I(is,h;ti) 
and I(v : h; £2) are compatible observables for any choice of the times; the same holds in 
the Heisenberg picture. Then, a theory of continuous observation of /(i/, h; t), t > 0, 
can be developed ]l2], |l4| and a probability law for the stochastic process "output 
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current" can be constructed; however, to obtain the spectrum we do not need the full 
theory, but only the second moment of /. Here and in the following for the expectation 
of any operator M we shall use the notation 

«M» = E [T W {M(p ® |e(/))(e(/)|)}] . (46) 
3.1.2. The power spectrum. In the long run the mean of the output power is given 

by 

P(v,h)= lim ^ [ T ((U(T)^I(u,h;t)) 2 U(T)))dt; (47) 

hi > has the dimensions of a resistance, it is independent of is, but it can depend 
on the other features of the detection apparatus. As a function of u, P{v, h) gives the 
power spectrum observed in the "channel h" ; in the case of the choice (|44|) it is the 
spectrum observed around the direction (0, 0) and with polarization e. 
The mean power P(v, h), under condition (fl5|), can be expressed as 



P(y,h) = , ^ + klk 2 Z(x;h), (48) 

47T 



fr r r~ r~ n 

X(x;h) = —J dte-«*Tr H ^D(h)U Kt [D{h) Poc \j +c.c, (49) 

v-^v 7 7 + 2/ , v x , Kn s 

D(h)=c- iS D(h)/(an), (51) 

|a| 2 /C[p] = £bM - BP +P P_ + BP_ P P+ ■ (52) 



D(h) is given by equation @, Poo by @, L B by Q-Q, y and z by @. 

In the expression jiq ) of the power the term fcjfc 2 /(47r) appears, independent off; 
it is a white noise contribution due to the detection scheme, coming out mathematically 
from the canonical commutation relations for the fields and known as "shot noise". 
The semigroup exp[/Ct] is trace preserving, but not positivity preserving, while p^ is 
the equilibrium state of Cb, which is a bona-fide Liouvillian, because it can be written 



in the Lindblad form ( g8c| ), ( ggjj ) 



3.1.3. Sketch of the proof of equations fclq) and (^|). The computations are similar 
to those of reference [[l9| , but now not only QSC but also classical stochastic calculus 
plays a role. Here we give only the main intermediate steps; a more detailed account 
can be found in p0| , where however the degeneracy of the levels is not considered. 

The first step is to rearrange in normal order the field operators involved in (f47|); 
in this way the constant shot noise appears and the spectrum reduces to 

h) = ^^((^ dAt ut (h; t) J* dA ont (h; s) e^^)^) )) + c.c. (53a) 



= dA lt(h; t)e- ivt jT dA out (h; s)e^ e~^)) > 0; (536) 
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the second expression shows the positivity of E(a;;/i), while the first one is useful 
for the successive developments. The output fields are given by A out (h;t) = 
U(t)iA(h;t)U(t) and we have 



dA out (h; t) = Y^(h\h k )U(ty fa S kidMt) + R k dt) U(t) 



(54) 



The second step is to show, by using (|54j), (|25|), fl2S|), (|45j) and the factorization 
properties of Fock space, that 



E(a;; h) = lim 



1 



x E 



dt 



ds exp < — 



+ \{v - uj) 



(t-s) 



cxp [v/B[W(t) - W(s)] } Tr H9 r{D(h)tC(t, s)D(h) 
(p(a)®|e(/)>(e(/)|)C(*,a)+}] + c.c. , 



C(t, s) = exp \ - ojt + yf§W{t) (P + - P_)\u(t)U(s)^ 



x exp 



VBW(s)] {P+ - P_) 



(55) 



(56) 



Finally an interplay between classical stochastic calculus and the quantum 
regression theorem, which holds for a dynamic like C(t, s), gives the expression ((fit). 



3.1. 4-. The angular distribution of the spectrum. Let us stress that is supported 
by the two states |1) = \F+,F+), |2) = \F-,F—) and that all the operations involved 
in (^) leave the span of |1), |2) invariant. So we are left with computations involving 
2x2 matrices; from now on we suppress the index Ti. in the trace and we use the Pauli 
matrices cr±, o~ Zi ... In particular we have P± = ^ (1 ± a z ) and 

D(h) = (h\l, 1; +1) £>i + (%+>P+ + (%->P- , (57) 

Di = - c~ i5 a- - \c lS + sin 6+ P+ - ie i5 ~ sin 5- P_ . (58) 

To obtain the angular dependence of the spectrum, let us introduce the states 
h± concentrated around (0, </>) and with right/left circular polarization, given by 
equation ( [51) with e — e±, where e±(6, <jrf — cx p^_"^) sm ^ zp cos q cqs,^ icos</>=F 

cos 9 sin 0, ± sin 0) . Then, the two angular spectra £± (x; #) , which, by the cylindrical 
symmetry of the problem, do not depend on 4>, are given by 

1 02 r+oo 

Z±(x;e) = —-X(x;h ± ) = —J^ dt C -^TT{D ± (6ye Kt [D ± (9) Poo }}+c.c., (59) 

D ± (6) = ± ± cos0) + J] <? e (0; ±)P . (60) 

The functions g e (9;±) — (h± \g e ) f yf\ AT| depend on 5% but after all they are free 
parameters of the theory: they are square integrable ^-functions, satisfying the 
constraint sin 6 [(1 + cos 6) ff e (fl; + ) — (! — cos 6) g e (Q; —)] d6 — 0, coming from the 



orthogonality of g e to |1, 1; +1), sec (306) 
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Then, by integrating over the whole solid angle, one gets the total spectrum 

E(x)=f d(9sin6> / d<p(Y,+ (x; 9) + E_ (x; 9)) = V T,(v; h k ) , (61) 
Jo Jo ' k 

where {hk} is any c.o.n.s. in Z. When one has g e = 0, as in the usual case, one gets 

3 /l±cos6» N 2 



Z±(x;6) = — { ) E(i). (62) 

For g e ^ the x and 9 dependencies do not factorize. In the experiments one measures 
something proportional to E+(a;; 9) + E_(x; 9) for 9 around 7r/2; this quantity fails to 
be proportional to E(x) only by the presence of some terms which we expect to be 
small and which are not qualitatively different from the other terms in E(x). So, for 
simplicity, we shall study only the total spectrum (|6l|). 

3.2. The total spectrum 

By choosing in ( |6l[ ) a basis with hi = |1,1;+1), h 2 = Ag/\\Ag\\, h 3 = 
\\Ag\\-i [||.g-|| 2 ||A g || 2 - | (Ag\g.) | 2 ] ~ 1/2 [||A. 9 ||V - {Ag\g-) Ag] , we get 

^Trj.D^e } = Tr \p\c Kt [D lPoo ]} + (g-\Ag) Tr{P +Poc } 

k 

+ (A.g| 5 _)Tr{P+e' c * [ Poo ]} + \\g-\\ 2 + \\ Ag\\ 2 Tr {P + e Kt [P +Poo ]} ■ (63) 
Then, from (fl9|), (|6l|), (|63|), we get the final expression of the total spectrum 



n 2 

Z7T 



- (w 3 + ie lS - sinS- + ifi 2 Wie ls sins) (d 2 - ie I(5 - sin<5_ 
q 

- ifi 2 die IS sin s) + e~ I<5 - sin 5 + (Sl 2 ci sin s - ie ls c 3 ) 

+ (d x + id 3 e ls sins) (u 3 + iQ 2 uie~ ls sins) + - {g~\g~ + £l 2 diAc/) 

n 2 

+ n 2 {Ag\g+} a - ft 2 d 3 || Ag\\ 2 u x + — {AgU + tfdiAg) v x + c.c. 

Q 





(64) 



v= KT- q w > c= kT^' u= KT- q [° 1 \> (65) 

/ 

K = G + y 1 | . (66) 
\ -1 

By integrating over the reduced frequency x, we get the strength of the spectrum 

-+oo 

E(x)dx = fi 2 |di (1 + Q 2 sm 2 S + ) + (1 - tt 2 di) (sin 2 <5_ + ||5_|| 2 ) 

+Re [da (e 2i5 - - l)] + n^Hs+f}. (67) 

Let us recall that |a| 2 is the natural line width, fl 2 is proportional to the laser 
intensity, z = (ui — uiq) /\a\ 2 is the reduced detuning, to which also the parameter 
z = z — fl 2 e is linked, y = B/\a\ 2 is the reduced laser bandwidth (|3lj|), O 2 1 a; | 2 e is an 
intensity dependent shift, x = [y — ojq) /\a\ 2 and 7 = 7/|a| 2 are the reduced frequency 



Quantum stochastic models of two-level atoms 



12 



and the reduced instrumental width, respectively, and q = i(x — z) + (7 + y) /2 (pO ) 



The other quantities entering all these equations are given by the equations (30£), 
( [30 <\ ), (|35|)-(|4l|); in particular, s = 5 + — 5_ and Ag = g + — g_. Let us note that 
£ i ll<?±l| 2 i (.9+ Iff-) are parameters linked to the scattering matrices, satisfying 
the constraints "|(<7+|<7_)| < \\g+\\ and "||Ag|| =0 =^> e = 0" (apart from 

this relation e is an independent parameter of the model). One can check that the 
spectrum T,(x) is invariant under the transformation: x — > —x, z — > —z, e — > —e, 
S± -> -6±, (g-\g+) -> (g+\g-)- 

3.2.1. The case of a monochromatic laser, y — 0. When the laser bandwidth vanishes 
the spectrum (p4 ) reduces to the one obtained in jy| , where the polarization of light 
and the level degeneracy were not considered; if also S ± = 1 is taken, it reduces to 
the Mollow spectrum jlq] . A particularity of the y = case is that the spectrum can 
be decomposed into an elastic and an inelastic part, which in turn are proportional to 
the elastic and inelastic electromagnetic cross sections, as explained in ||l9| ; for y > 
such a decomposition has no meaning. Figures 2 and 3 of E| compare the spectra 
(for 7 = 0.6 and for various detunings and laser intensities) in the Mollow case and 
in the modified case with 6+ = -0.03, 5- = 0.13, ||.g±|| 2 = 0.005, (g+\g-) = -0.005, 
e = —0.001; the two figures have x — z in the horizontal axis and E/f2 2 in the vertical 
one. The parameters have been chosen in such a way that the modified spectrum 
would not too different from the Mollow one in the resonant case. In plotting the 
spectrum a factor 1/2 has been inserted by error in the inelastic part, so the two 
figures have only a qualitative meaning; in any case, one sees the symmetry of the 
Mollow spectrum which is lost in the modified case. 

3.2.2. The case = 1. Let us recall that the usual model, with only the 
absorption/emission process, corresponds to 5± = 0, g± = 0, e = 0, z = 2, s = 0. In 
this case we obtain 



Zn 



1 , /„ 40 



-v 3 d + \^2v 3 + — ) Red 



c.c. , (68) 



1 + y + 2iz o , , , 

d= 4 J +T2 , F 2 = (l + y)(l + y + 20 2 ), (69) 

v 3 = [2 + 7 + y + 2i (a; - z)] [1 + 7 + % + 2i (a: - 2z)] /N , (70) 

N = 40 2 [1 + 7 + 2y + 2i(x - z)] + [2 + 7 + y + 2i (x - z)] 

x [1 + 7 + Ay + 2i (x - 2z)] (1 + 7 + 2\x) . (71) 

Now, the spectrum T,(x) is invariant under the transformation: x — * — x, z — > — z. 

If we put also 7 = 0, which means that the instrumental width is negligible, then 
one can check that the fluorescence spectrum T,(x) coincides exactly (apart from the 
different normalization) with the spectrum computed by Kimble and Mandel |l(| [TtJ ; 
when also the laser bandwidth y vanishes, the equations above reduce once again to 
the Mollow spectrum |lj| . One of the more interesting results of Kimble and Mandel 
is that a not vanishing laser bandwidth renders asymmetric the Mollow spectrum and 
enhances the side peak near u>o, as shown in the figures of []l6| . In our notations, the 
figures 2, 3, 4 of |l(| give, for 7 = 0, irY, on the z-axis versus 2(x — z) on the x-axis 
and y on the y-axis; figure 2 is for fi 2 = 0.0025 and z = —2.5, figure 3 for £! 2 = 25 
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Figure 1. The spectrum for 7 = 0.6, C 2 = 28, 2 = and laser bandwidths 

» = , 0.5, 1, 4. Solid lines: modified model with the parameters given in section 
3.2.5; dashed lines: usual model. 



and z = 0, figure 4 for £1 2 = 25 and z = —1.5. The figures 5 and 6 of |ll| give, for 
7 = 0, 7rS on the z-axis versus 2(x — z) on the x-axis and 1z on the y-axis; figure 5 is 
for n 2 = 25,y = 0.5 and figure 6 is for ft 2 = 2500, y = 10. 



3.2.3. Low intensity laser. For a laser of low intensity the spectrum reduces to 

TtSOe) , si/ m2 A, ,,2 7 1 / m2\ 

2 ^ r22 =0 V 7 + y / 

+ y|c(x -2,7 + 2/)a(y, z) + c(-x, l+j)c(-z, 1 + y)| 2 
+ 7J/ |c(-x, 1 + 7 )| 2 \c(-z, 1 + j/)| 2 , (72) 

a(y, z) = c(-z, 1 + y) - ie"' 5 - sin S- , c(x, A) = A ^ 2ix ■ (73) 

In this case one sees well why a decomposition in elastic and inelastic parts is not 
possible for y > 0. Indeed, one has two peaks, one centered at v = u> (laser frequency) 
and one centered at v = luq (atomic frequency); the first peak can be interpreted 
as coming from an elastic channel and the second one as coming from an inelastic 
channel, but the second term in the r.h.s. of (|72|) gives rise to an interference term 
between the two channels. 
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Figure 2. The spectrum for 7 = 0.6, f2 2 = 28, z = 2.5 and laser bandwidths 
» = D , 0.5, 1, 4. Solid lines: modified model with the parameters given in section 
3.2.5; dashed lines: usual model. 



3.2.4- The broadband case. For a very large bandwidth, j/>i, 2, f2 2 , the spectrum 
becomes independent from the detuning and reduces to 

lim 7^ = \ c ( x ~VA + n)+ie [S - sinJ I 2 + k \c(x - rj, 1 + n)\ 2 + \\g_\\ 2 , (74) 

where c(a;, A) is given in ( |73| ) and 

77 = SI 2 fe- J sin2s^ , k = 7 + ft 2 (sin 2 s + || A.g|| 2 ) . (75) 

Note the 17-dependent position of the peak, which can be interpreted as a light shift. 
In the usual case S(x) becomes purely Lorentzian, centered in zero and with f2- 
independent width: 

Try I + 7 

^ ^ = (l+ 7 ) 2 +4x 2 • (76) 

3.2.5. Plots. In the general case the total spectrum is given by equations (piQ-diq). 
In figures 1, 2 and 3 we compare the spectrum predicted by the usual model (dashed 
lines), in which only the absorption/emission channel is present, with the spectrum 
predicted by our modified model (solid lines), in which both the absorption/emission 
channel and the direct scattering channel arc present. In both models we take 7 = 0.6, 
O 2 = 28; the usual model is characterized by S± — 0, ||g±|| 2 = 0, e = 0, while as an 
example of modified model we choose 5+ = —0.03, 6- = 0.13, ||.<?+|| 2 = 0.0045, 
||S_|| 2 = 0.0055, (g+\g-) = -0.004 + i x 0.002, e = -0.001. We consider three laser 
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Figure 3. The spectrum for 7 = 

bandwidths y = 0, .5, 1, 4. Solid lines: 



0.6, = 28, 2 = -2.5 and laser 
modified model with the parameters 



given in section 3.2.5; dashed lines: usual model. 



detunings z = 0, 2.5, —2.5 and four laser bandwidths y — 0, 0.5, 1, 4. Our choice 
of the parameters for the modified model is such that for a monochromatic laser in 
resonance (figure 1(a)) the modified spectrum is not quantitatively too different from 
the usual one, but its asymmetry is clear. The differences between the two cases are 
enhanced by the presence of both the detuning and the bandwidth. 

Finally in figure 4 we consider a large bandwidth; ho wever , y is not much larger 
than SI 2 , so that we are not exactly in the case of section 3.2.4 . Anyway, only a peak 
survives and the spectrum is nearly independent from the detuning. In the usual case 
we see a nearly Lorentzian shape, centered at x — (y near the atomic frequency u)q) 
and with a very weak residual ^-dependence. In the modified case the spectrum is 
strongly asymmetric and the position of the peak and the width are O-dependent. 
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